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Abstract 

We provide evidence for a high precision model-independent low momentum nucleon- 
nucleon interaction. Performing a momentum-space renormalization group decima- 
tion, we find that the effective interactions constructed from various high precision 
nucleon-nucleon interaction models, such as the Paris, Bonn, Nijmegen, Argonne, 
CD Bonn and Idaho potentials, are identical. This model-independent low momen- 
tum interaction, called Vjowk) reproduces the same phase shifts and deuteron pole 
as the input potential models, without ambiguous assumptions on the high mo- 
mentum components, which are not constrained by low energy data and lead to 
model-dependent results in many-body applications. V\o^ k is energy-independent 
and does not necessitate the calculation of the Brueckner G matrix. 

PACS: 21.30.Cb; 21.60.-n; 21.30.Fe; ll.lO.Hi 
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1 Introduction 



In low energy nuclear systems such as finite nuclei and nuclear matter, one 
can base complicated many-body calculations on a simple picture of point-like 
nucleons interacting by means of a two-body potential and three-body forces 
when needed. Unlike for electronic systems where the low energy Coulomb 
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force is unambiguously determined from Quantum Electrodynamics, there is 
much ambiguity in nuclear physics owing to the non-perturbative nature of 
Quantum Chromodynamics (QCD) at low energy scales. Consequently, there 
are a number of high precision, phenomenological meson exchange models of 
the two-nucleon force Vnn, such as the Paris [1], Bonn [2], Nijmegen [3], Ar- 
gonne [4] and CD Bonn [5,6] potentials, as well as model-independent but less 
accurate treatments based on chiral Effective Field Theory (EFT) [7,8,9,10], 
for a review see [11]. We also study the Idaho potential [12], which is based 
on the EFT framework, but some model dependence is introduced in order to 
achieve similarly high precision as compared to the other interaction models. 

These nuclear force models incorporate the same one-pion exchange interac- 
tion (OPE) at long distance, but differ in their treatment of the intermediate 
and short-range parts (e.g., parameterization of the repulsive core compared 
to uo exchange, different form factors, dispersive or field theoretical treatment 
of the 2tt exchange). The fact that the different short distance constructions 
reproduce the same phase shifts and deuteron properties indicates that low 
energy observables are insensitive to the details of the short distance dynam- 
ics. The EFT approach exploits this insensitivity by explicitly keeping only 
the pion and nucleon degrees of freedom (in accordance with the spontaneous 
breaking of chiral symmetry in the QCD vacuum) and encoding the effects of 
the integrated heavy degrees of freedom in the form of couplings which mul- 
tiply model-independent delta functions and their derivatives, see e.g. [13]. 
The EFT thus provides a model-independent description of the two-nucleon 
system. However, the high precision, i.e., x^/datum k, 1, description of the 
nucleon-nucleon scattering data provided by conventional models is at present 
not achieved by the rigorous EFT potentials. 
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In Fig. 1, we observe that the rcahstic models of Vnn have quite different 
momentum-space matrix elements despite their common OPE parts and re- 
production of the same low energy data. It demonstrates that the low energy 
phase shifts and deuteron properties cannot distinguish between the models 
used for the short distance parts. However, in many-body calculations the 
assumed short-distance structure of a particular Vnn enters by means of vir- 
tual nucleon states extending to high momentum, which will lead to model- 
dependent results, e.g., in Brueckner Hartree-Fock calculations of the binding 
energy of nuclear matter. It is clearly of great interest to remove such model 
dependence from microscopic nuclear many-body calculations. 

Motivated by these observations, we perform a renormalization group (RG) 
decimation and integrate out the ambiguous high momentum components 
of the realistic interactions. Due to the separation of scales in the nuclear 
problem, it is reasonable to expect that the model dependence of the input 
potentials will be largely removed as the high momentum components are 
excluded from the Hilbert space. Starting from any of the Vnn in Fig- 1, 
we integrate out the momenta above a cutoff A to obtain an effective low 
momentum potential, called V\o^ k- The physical condition is that the effective 
theory reproduces the deuteron pole and the half-on-shell (HOS) T matrices 
of the input Vnn rnodel (i.e., the observable scattering phase shifts and the 
low momentum components of the two-body wave functions, which probe the 
OPE part of the interaction), but with all loop momenta cut off at A. We 
will find the striking result of Fig. 2 that the V^iow k becomes independent of 
the particular input model for A < 2.1 fm~^, which corresponds to laboratory 
energies Ei^h < 350 MeV. The latter is the energy scale over which the high 
precision potentials are constrained by experimental data. 

Our strategy thus is a compromise between conventional models and EFT 
treatments. The resulting Viow k reproduces the nucleon-nucleon scattering 
phase shifts with similar accuracy as the high precision potentials, but without 
making further assumptions on the detailed high momentum structure, which 
cannot be resolved by fitting to the low energy data only. 



2 Renormalization group decimation 



The first step in the RG decimation is to define an energy-dependent effective 
potential (called the Q box or the Bloch-Horowitz potential in effective inter- 
action theory), which is irreducible with respect to cutting intermediate low 
momentum propagators. The Q box resums the effects of the high momentum 
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modes, 



Q[k ,k]u) ^V^^[k ,k) + -V p dp, (1) 

TT JA UJ — P'^ 



where k' , k, and p denote the relative momentum of the outgoing, incoming, 
and intermediate nucleons. The principal value HOS T matrix for a given 
partial wave is obtained by solving the Lippmann-Schwinger equation 

^fW u u2\ T/ fu u\^'^ r ^NN(A:^p) T{p,k;k'^) ^ 

T{k ,k]k ) = V^^{k ,k) + -V 5 pap. (2) 

t: Jo k'^ — 

In terms of the effective cutoff-dependent Q potential, the scattering equation 
can be expressed as 

T(fc', k- k') = Q{k', k- e) + -V r ^^^'^P'^f^ TiP,k;e) 

n Jo k'^ — p'^ 

The low momentum effective theory defined by Eq. (1) and Eq. (3) preserves 
the low energy scattering amplitudes and bound states independently of the 
chosen model space, i.e., the value of the cutoff A. However, the energy de- 
pendence of the effective Q potential is inconvenient for practical calcula- 
tions. In order to eliminate the energy dependence, we introduce the so-called 
Kuo-Lee-Ratcliff folded diagrams, which provide a way of reorganizing the 
Lippmann-Schwinger equation, Eq. (3), such that the energy dependence of 
the Q{k' , k;uj) box is converted to a purely momentum dependent interaction, 
Mow k{k', k) [14,15] . The folded diagrams are correction terms one must add to 
Eq. (3), if one were to set all Q box energies right side on-shell. This explicitly 
leads to (for details see [15,16]) 



Mow k(A;', k) = Q{k', k-e)+ - V f p^dp Q{p, k; P 

TT Jo 

Q{k',p; e)-Q ik',p;p') 
k"^ — pi 



The folded diagram resummation indicated in Eq. (4) can be carried out to all 
orders using the similarity transformation method of Lee and Suzuki [17,18]. 
By construction, the resulting energy-independent Mow k preserves the low 
momentum HOS T matrix of the input Vnn model [19] , 

T{k', k; e) = Mow k) + -vt ^-^'^^'^ Tij^JrJ,^) 

t: Jo k^ — p^ 

where all momenta are constrained to lie below the cutoff A. As our RG 
decimation preserves the half-on-shell T matrix, we have dT{k', k; k'^)/dA = 0, 
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Fig. 2. Diagonal and off-diagonal momentum-space matrix elements of the Vjow k ob- 
tained from the different bare potentials in the ^Sq channel for a cutoff A = 2.1 fm"^^. 

which implies a RG equation for Viowk [16] 



(6) 



Similarly, a scaling equation is obtained for the Q box by integrating out an 
infinitesimal momentum shell, and one finds 



_d_ 

dA 



2Q{k',A;p^)Q{A,k;p^ 



TT 



(7) 



This equation was obtained previously by Birse et al. by requiring the invari- 
ance of the full-off-shell T matrix, dT{k\ k;p^)/dA = [20]. 

The RG equation, Eq. (6), lies at the heart of the approach presented here: 
given a microscopic input model with a large cutoff, one can use the RG 
equation to evolve the bare interaction to a physically equivalent, but simpler 
effective theory valid for energies below the cutoff. The RG evolution separates 
the details of the assumed short distance dynamics, while incorporating their 
detail-independent effects on low energy phenomena through the running of 
the effective interaction [16,21]. 

The principal difference between the presented RG approach and the standard 
EFT one is that we do not expand the interaction in powers of local opera- 
tors. This implies that, first, one does not make assumptions on locality, and 
that second, we do not truncate after a certain power included in the low mo- 
mentum interaction. In the EFT approach, power-counting is used to calculate 
observables to a given order with controllable errors. In contrast, we start from 
a Hamiltonian in a large space, which reproduces the low-energy observables 
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with high accuracy and then decimate to a smaller, low momentum space so 
that the observables are reproduced. Since the exact RG equation is solved 
without truncation, the invariance of the T matrices is guaranteed. 

Nevertheless, there are close similarities between the presented RG approach 
and the standard EFT one. For example, at cutoffs above the pion mass, V\ow k 
keeps the pion exchange explicit, while the unresolved short distance physics 
could be encoded in a series of contact terms. Moreover, we find below that 
V\ow k seems only to depend on the fact that all potential models have the same 
pion tail, and fit the same phase shifts. This is similar to EFT treatments, 
where the interaction is constrained by pion exchange, phase shifts and the 
choice of regulator. 



3 Results 

Referring to Fig. 2, we find the central result of this Letter. The diagonal ma- 
trix elements of Viow k obtained from the different Vnn of Fig. 1 collapse onto 
the same curve for A < 2.1 fm~^. Similar results are found in all partial waves 
and will be reported elsewhere in detail [23] . We emphasize that Viow k repro- 
duces the experimental phase shift data and the deuteron pole with similar 
accuracy as compared with the high precision potential models, as we show 
in Fig. 3. The reproduction of the phase shifts with the renormalized V^iow k 
does however not require ambiguous high momentum components, which are 
not constrained by the low energy scattering data. We note that the change 
in sign of the ^Sq and the ^Pq phase shifts indicates that the effects of the 
repulsive core are properly encoded in Viow k- 

Intriguingly, the Viow k mainly differs from the bare potential by a constant 
shift. This was previously observed by Epelbaum et al. for a toy two- Yukawa 
bare potential [24]. The constant shift in momentum space corresponds to a 
smeared delta function in coordinate space and accounts for the renormaliza- 
tion of the repulsive core from the bare interactions, see also [25]. 

In Fig. 2, we find a similar behaviour for the off-diagonal matrix elements, 
although the T^ow k derived from the Idaho potential begins to differ from the 
others at approximately 2 771,^ = 1.4 fm^^. This discrepancy in the off-diagonal 
matrix elements arises from the fact that the Idaho potential used here treats 
the 27r exchange differently than the meson models do. We can integrate out 
further and lower the cutoff to A < 1.4 fm"^, then the off-diagonal elements 
collapse as well. 

These results can be understood from T matrix preservation. The HQS T 
matrix determines the phase shifts as well as the low momentum components 
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Fig. 3. S-wave (singlet and triplet with mixing parameter) and P-wave phase shifts of 
^low k for a cutoff A = 2.1 fm~^ compared to the input Vnn model used, here the CD 
Bonn potential. We also show the results of the multi-energy phase shift analysis 
(PWA93) of the Nijmegen group [22]. For Viow k the agreement of the calculated 
phase shifts with the PWA93 is determined by the quality of the fit for the bare 
interaction. 

of the low energy scattering and bound state wave functions. Using the spectral 
representation of the T matrix, Viow k can be expressed as 

V^ow k(A;', k) = T{k', k; k^) + - V r T{k', p; T{p, k; (8) 

TT Jo — /e^ 

The Vnn models give the same on-shell T matrices over the phase shift equiv- 
alent kinematic range of < 2.1 fm~^, but their off-shell behavior is a priori 
unconstrained. In practice, however, one observes that the realistic potential 
models result in similar HOS T matrices at low energies and momenta. This 
is understood from the fact that these models share the same long-range OPE 
interaction and differ most noticeably on short- distance scales set by the repul- 
sive core, r ~ 0.5 fm and smaller. Consequently, one expects that the off-shell 
differences are suppressed at energies and momenta below a corresponding 
scale of A ~ 1/r ~ 2.0 fm~^. It is clear from Eq. (8) that the approximate 
HOS T matrix equivalence is a sufficient condition for the Viow k to be inde- 
pendent of the various potential models. Moreover, the spectral representation 
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Fig. 4. Evolution of the diagonal momentum-space matrix elements for the Vjow k 
derived from the different bare potentials from larger values of the cutoff to 
A = 2.1 fm~^ in the ^Si partial wave. 

clarifies that the off-diagonal matrix elements of V\o^ k are more sensitive to 
a particular off-shell behavior as observed in Fig. 2, and thus deviate at a 
lower cutoff than the diagonal matrix elements. The collapse of the diagonal 
momentum-space matrix elements at the scale set by the constraining scat- 
tering data, A ^ 2.1 fm~^, is nicely illustrated in Fig. 4, which shows the RG 
evolution as the cutoff is successively lowered in the '^Si partial wave. 



Next, we analyze the scaling properties of l^iow k- For this purpose, we show the 
matrix element Viowk(0,0) versus cutoff in Fig. 5 for the ^Sq and ^Si partial 
waves. The main results are the following. T matrix preservation guarantees 
that V\o^ k(0, 0) flows toward the scattering length as the cutoff is taken to 
zero, these are aiso = —23.73 fm and aasj = 5.42 fm and are well reproduced 
with V"iow k- 

For cutoffs A > rriT,, l^owk is nearly independent of the cutoff in the ^Sq 
channel and weakly linearly dependent in the ^Si channel. According to EFT 
principles, the couplings are nearly independent of the cutoff as long as A 
is large enough to explicitly include the relevant degrees of freedom needed 
to describe the scale one is probing. Thus, the running of l^ow k initiated at 
A ~ m^r is a result of integrating out the pion. The rapid changes at very 
small A are the result of the large scattering length. To illustrate the effects of 
the large scattering length, we solve the RG equation for V^ow k(0, 0) for small 
cutoffs, with the ^Sq scattering length as the A = boundary condition. One 
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Fig. 5. RG flow of yiowk(0,0; A) versus cutoff A in the ^Sq and ^Si partial waves. 
The solid line represents the solution of the RG equation for small A as discussed 
in the text. 



finds 

^iowk(0,0)^-— ^— -, (9) 
1/aiSo - 2A/7T 

which is shown as the solid curve in Fig. 5. Clearly, the agreement for small 
A < m7r/2 is convincing. Eq. (9) is identical to the Kaplan-Savage- Wise so- 
lution for the renormalization of the momentum-independent contact term in 
the pionless EFT, upon identifying A with the regulator mass fi used in di- 
mensional regularization [8]. The weak dependence on the cutoff in the ^Si 
partial wave results from the dominantly second order tensor contributions, 
which are peaked at relative momenta k ^ 2 fm~^ [26]. 

Finally, we note that V^ow k preserves all scattering and bound states with 
energy \E\ < A^. Therefore, the jump in the ^Si channel is the bound state 



contribution to the T matrix and occurs at A = /cq = ^JmED/K^ = 0.23 fm~^. 
From effective range theory and the Low equation [27], we have for the bound 
state contribution and thus the discontinuity 



Al^ow k(0, 0; M3Si = ; ^ = 22.8 fm, (10) 



where ro = 1.75 fm denotes the effective range in ^Si. This is in very good 
agreement with our results. 
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4 Summciry 



In this work, we have shown that the separation of scales in the nuclear force 
can be successfully applied to derive a model-independent low momentum 
nucleon-nucleon interaction. Viowk is obtained by integrating out the high 
momentum components of various potential models, leading to a physically 
equivalent effective theory in the low momentum Hilbert space. The decima- 
tion filters the details of the assumed short-distance dynamics of the bare 
interactions, provided one requires that the low energy observables are pre- 
served under the RG. Wc have argued and provided numerical evidence that 
the model-independence of Viow k is a consequence of the common long-range 
OPE interaction and the reproduction of the same elastic scattering phase 
shifts. The momentum scale A ~ 2.1 fm~^, corresponding to laboratory en- 
ergies E'lab ~ 350 MeV, is precisely the scale at which the low momentum 
interaction becomes independent of the input models. Our work demonstrates 
that the differences in the high momentum components of the nucleon force 
are not constrained by fits to the low energy phase shifts and the deuteron 
properties. 

Mow k does not have high momentum modes, which are related to the strong, 
short range-repulsion in conventional models. Therefore, the low momentum 
interaction is considerably softer than the bare interactions (both in a plane- 
wave and an oscillator basis) and does not have a repulsive core. As a conse- 
quence, when Viow k is used as the microscopic input in the many-body prob- 
lem, the high momentum effects in the particle-particle channel do not have to 
be addressed by performing a Brueckner ladder resummation or short-range 
correlation methods. In fact, with a cutoff on relative momenta, the phase 
space in the particle-particle channel is comparable to the particle-hole chan- 
nels, and it would seem strange to resum the particle-particle channel, while 
the particle-hole channels are treated perturbatively in a hole-line expansion. 

The use of Viow k in microscopic nuclear many-body calculations leads to 
model-independent results. VJow k has been successfully used as shell model 
effective interaction in model space calculations for two valence particle nuclei 
such as ^^O and ^^^Te [28]. The starting point of these calculations has tradi- 
tionally been the Brueckner G matrix, which depends on the bare Vnn used 
as well as the particular nuclei via the Pauli blocking operator. By means of 
Viowk, the same low momentum interaction is used in different mass regions. 
Viowk has further been incorporated into Fermi liquid theory. This connects 
the low momentum interaction in free space and the quasiparticle interaction 
in normal Fermi systems. Two constraints have been derived which relate the 
Fermi liquid parameters of nuclear matter to the S-wave low momentum inter- 
action at zero relative momentum [29]. Finally, adopting the RG approach to 
Fermi systems proposed by Shankar, Vjow k is taken as the input to microscopic 
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calculations of the quasiparticle interactions and the pairing gaps in neutron 
matter [30] . From the success of these applications, we believe that the model- 
independent V^ow k is a very promising starting point for microscopic nuclear 
many-body calculations. 

Few-body calculations with V\ow k are in progress. Preliminary results for the 
ground state energies of triton and ^He show that the three-body force for 
Viowk is weaker than e.g., the conventional three-body force constructed for 
the Argonne potential. 
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